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Abstract. The Hannay angle has been previously studied for a celestial circular 
restricted three-body system by means of an adiabatic approach. In the present 
work, three main results are obtained. Firstly, a formal connection between 
perturbation theory and the Hamiltonian adiabatic approach shows that both 
lead to the Hannay angle; it is thus emphasised that this effect is already 
contained in classical celestial mechanics, although not yet defined nor evaluated 
separately. Secondly, a more general expression of the Hannay angle, valid for 
an action-dependent potential is given; such a generalised expression takes into 
account that the restricted three-body problem is a time-dependent, two degrees 
of freedom problem even when restricted to the circular motion of the test body. 
Consequently, (some of) the eccentricity terms cannot be neglected a priori. 
Thirdly, we present a new numerical estimate for the Earth adiabatically driven by 
Jupiter. We also point out errors in a previous derivation of the Hannay angle for 
the circular restricted three-body problem, with an action-independent potential. 
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1. Introduction 

During the 20 years since Berry £Q discovered the quantum geometric phase and 
Hannay [2] demonstrated its classical counterpart, a great deal of research effort has 
gone into the subject of geometric phases. Shortly after their discovery, Berry [3] 
showed how they were formally related. The Hannay effect is known to occur in a 
system undergoing an adiabatic change. The action variables behave as adiabatic 
invariants, which remain unchanged as the system evolves. The conjugate angles, 
however, being affected by the system's slow variation, exhibit a geometric shift which 
arises from the Hamiltonian being cyclic. 

From the mid 80s onwards, studies devoted to the Hannay angle, have analysed several 
applications, including measurements in the solar system related to the non-sphericity 
and slow rotation of the Earth 0]. Instead, Berry and Morgan |Hj have attempted 
to determine the Hannay angle for the restricted circular three-body problem. They 
give estimates for the Earth and Mars revolving around the Sun, under the adiabatic 
influence of Jupiter, and for a geostationary satellite, under the influence of the Moon. 
We have structured the paper in the following manner. In the second section, we 
derive the Hannay angle for the restricted circular three-body problem based on 
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classical perturbation analysis of celestial mechanics, using the same simplified action- 
independent potential Hamiltonian of [5]. This derivation simply underlines that the 
Hannay effect is not a new phenomenon in celestial mechanics, but has always been 
taken into account in all high order analytic theories. However, never before has it 
been defined nor evaluated separately. Furthermore, the previous derivation of the 
Hannay angle [5] contains an error, (with reference to the masses of a three-body 
system, we let the attractor be A, the perturber P and the test mass be T) due to 
forcing T to move on a fixed, circular orbit. As a result, the interaction between P 
and A is neglected. In the third section, we start from a proper representation of 
the restricted circular three-body problem, with an action-dependent potential and 
derive a generalised Hannay angle. We show that, although we are concerned by 
T circular motion solely, some eccentricity terms cannot be neglected a priori, as the 
restricted three-body problem is intimately a time-dependent problem characterised by 
two degrees of freedom. Neglecting all eccentricity terms as in [S] , implies an artificial 
reduction to one degree of freedom (as if T were forced by an external constraint to 
move on a circular guide). Finally, the appendix gives some useful basic expressions 
of the three-body problem. 



2. The simplified restricted circular three-body problem 

2.1. Adiabatic and perturbative approaches 

Berry and Morgan have determined the Hannay angle for a simple system, inspired 
by the restricted circular three-body problem, with the following Hamiltonian: 

H(q,p;t) = ^ P 2 + V cos(q~et) , (1) 

where Vq ~ Mp/Ma is a small parameter, of the order of the mass of P relative to A, 
and e is another small parameter, proportional to the ratio of the orbital frequencies 
of P and T. The Hamiltonian is not the Hamiltonian of the restricted three-body 
problem, which we shall indeed consider in the next section. 

Within the framework of the problem described by eq. and by the adiabatic theory, 
the Hannay angle was approximated, at lowest order in Vq by [5]: 

where /, the action, is the adiabatic invariant: 

I=±- Ipdq = J- / y/2(E-V cos q)dq (3) 
2ir J 2tt J q 

E is the T kinetic energy and v = dE/dl is the angular frequency. The Hannay 
effect is defined if closed curves of constant action return to the same curves in phase 
space after a time evolution. It is the extra angle picked up by the angle variables, in 
addition to the time integral of the instantaneous frequencies. We remark that, given 
that the expression of the Hannay angle in (2) is negative, it causes a lengthening of 
the orbital period, conversely to what is stated above eq.(49) in [II 
We now re-derive the expression of the Hannay angle J2J for the Hamiltonian JJJ but 
using the approach of classical perturbation theory. The first step is to re-write the 
Hamiltonian Q in autonomous form. To this purpose, we introduce a pair of conjugate 
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s x f = / + {/,*} + M,xhx} + ••• , (6) 



action angle variables T ,t = and, in order to have f = e in the Hamiltonian 
equations, write as: 

H{q, p, r, T) = ip 2 + eT + V cos(g - r) . (4) 

We now consider ip 2 + eT as an unperturbed, integrable Hamiltonian Hq and 
Vo cos(q — r) as a small perturbation Hi, of order 1 in Vq. Classical perturbation 
theory consists in finding a change of variables (q,p,r,T) — > (q' ,p' ,t' ,T') so as to 
remove the term 7ii from the new Hamiltonian (i.e. from the Hamiltonian expressed 
in the new variables). The simplest way to achieve this result is to introduce the 
new variables by Lie series (see ch. 2 of jSJ). Hence, we look for a transformation of 
variables of the form 

P = S x p' , q = S x q' , T = S X T' , r = S x t' , (5) 

where S x is the Lie series operator, which, applied to a generic function /, reads 

1 

T 

where {., .} denotes the Poisson bracket and x is an a priori unknown function. The 
introduction of the new variables through (JjjJ has the advantage that (i) the new 
variables are canonical by construction, and (ii) the Hamiltonian can be expressed in 
the new variables in a straightforward manner, as 

1 

2 

From the above expression, for the elimination of the term of order 1 in Vq from TC , 
the function x must satisfy the equation 

{Ho,x} + Hi = Q. (8) 

In our case, given the expression of Tlo and Hi, we need to set 

X = r^—sin(q - t) . (9) 

p — e 

Now that the function \ has been defined, terms of H 1 of order 2 in Vq can be computed 
from Q) and are: 

\{{n Q , x}, x} + {ni, x } = 2 J[ e) 2 sin V ~ O ■ ( 10 ) 

So, the new Hamiltonian up to order 2 in Vq reads: 

n ' = l p ' 2+er+ w^ s[n2{q '- T,) - 

We now average Ijllf) with respect to q' . The averaged Hamiltonian reads 
The equation of motion for q' is 



H' = S X H = H + {H, X } + xh X} + ■ ■ ■ (7) 
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Referring to the slow periodic variation in time as modulation, this derivative can be 
split into an un-modulated frequency (evaluated for e = 0), namely p' — V 2 /[2(p') 3 ], 
and a correction, obtained by derivation of (|13|l 

Sq' = -3eV 2 /[2( P r] ■ (14) 

The Hannay angle is defined as the shift of q' over a modulation period 27r/e with 
respect to the un-modulated case, and is therefore 



2 



»— t"— w - - (15) 

The expressions |PJ and (|15|l are formally identical, but the former is a function 
of I while the latter is a function of p' . We now show that / and p 1 coincide, up to 
order O(V 2 ) and 0(e). In fact, expanding © as a Taylor series of Vo and recalling 
from JU that E = p 2 /2 + Vo cos(q — et), we obtain: 

I = p+— cos{q - et) + O{V 2 ) . (16) 
V 

On the other hand, applying the transformation JSJ) with the generating function \ 
given in JJJJ, we get: 

p' = P + cos(q - r) + O(V 2 ) . (17) 

Therefore p' = I up to order 2 in Vo and order 1 in e. This result, together with (JSJ 
and l|15|l proves that the adiabatic approach and the perturbation approach lead to 
an identical Hannay angle, at least for the Hamiltonian used in 



2.2. Hannay angle in the simplified problem 

The derivation in {jy forces the perturbed body to move on a circular orbit around 
the Sun, rather than considering from the beginning the restricted three-body problem 
(RTBP). In this section we compute the Hamiltonian that substitutes (1), if the RTBP 
is considered without constraints. 

As shown in the Appendix, the non Keplerian part of the potential in a three-body 
problem (to be added to the Keplerian AT term) is given by||: 

^ G K«^t;" s2 Ht;) <18) 

where R is the distance between bodies and a is the angle centred on A between T 
and P. We note that there is a cancellation of the linear terms in cos a, that was not 
previously considered, eq.(A4) of see appendix^ 

|| The potential in |5] has the dimensions of a frequency squared. Herein we use conventional units. 

The source of the error in |5] is caused by forcing the Earth to move on a fixed curve (a circle), 
with two consequences. The former is the absence of the mutual attraction of P and A, second term 
of ea. 1551 . and thus the non-cancellation of the first order term due to the action of P on T, eg. 1561 . 

As a result, the total action of P is proportional to ^ ( fl^p ) ms tead of R ^ p \ r TA ) as m ec l- 
(43) of 0. The latter is to disregard a priori the potential terms proportional to the eccentricity, i.e. 
the removal of a degree of freedom. These terms induce a modification of the Hannay angle of order 
in the eccentricity, eq. 1421 and do remain even if we set e=0 a posteriori. 
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Since the term l/R^p does not depend on T coordinates, the useful part of the 
perturbative function reduces to: 

where we have used a = \t — \p, At and Xp being the true longitudes of T and P 
respectively. The orbit of P is supposed to be circular; therefore Rap is constant and 
the true longitude Xp is equal to the mean longitude. But, as proved in section 3, for 
the orbit of T the circular approximation is not valid. The perturbation needs to be 
expanded up to the first power of the eccentricity e of T. Also, for an eccentric orbit 
the true longitude and the mean longitude are different, and it is more convenient to 
write the perturbation as a function of the latter. Using the classical expansions of 
the elliptic motion up to first power of the eccentricity 



Rta 



1 - e cos(A Tm + tp) (20) 



we get: 



cos(2At) = cos(2ATm) + 2e [cos(3ATm + tp) — cos(ATm — tp)] (21) 
sin(2A T ) = sin(2A T m) + 2e [sin(3AT m + tp) - sin(A T m - tp)} (22) 

y , = GMpf a \ l {1 + 3cos2(ATm _ Ap)+ (23) 
Rap V Map J 4 

e [—2 cos(ATm + tp) — 9 cos(A-rm — tjj — 2Ap) + 3 cos(3AT m + tp — 2Ap)]} 

where a is the semi-major axis of the orbit of T around A, Xxm is the mean longitude 
of T and tp is the opposite of the argument of its pericentre. 
We now introduce the actions 



A T = ^GM A a = At (l — y/l — e 2 ) (24) 



conjugate to the angle Arm an d tp respectively. We make use of the third Keplerian 
law relating the orbital frequency lot of T to the distance from A, 

2 GMa 

lo t = — — (25) 

Similarly we introduce the notation 
GM, 



and get 



L AP 



2 M P A T 1 



V = ^-f- - {(1 + 3 cos [2(A Tm - Xp)] + (27) 

IVLa w T 4 

-2cos(ATm + tp) — 9cos(ATm — tp — 2Xp) + 3cos(3Atto + tp — 2Ap)]| 

Now, V' depends on the actions and the formalism used by Berry and Morgan is no 
longer directly applicable. Instead, the Hannay angle will be derived by means of a 
perturbation approach in the next section. 

The expression (|19H of the perturbation is derived from an expansion of l/R A p in 
powers of Rta/ Rap limited to the second order. This expansion by means of Legendre 
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polynomials is inaccurate when Rta/ Rap is not small. In this Fourier like 

expansion, using Laplace coefficients 0, is more appropriate: 



Rta a 



Rap Rap 




1 - l/-T— COs(A T m + 1p) 

At 



a 



fe=i 




where the Laplace coefficients bf^ 2 are functions of 



potential takes the form: 



V 



n _ M P 



M A 



B k {A T ) cos[fc(A Tm - A P ) 



(28) 



COs(A T m + 1p) 



Rap GM A R A p 



A 2 

' Finally, the 



(29) 



C+(A T )V^cos[/cA rm - i> - (k + 1)A P ] + C fe -(A T )V*cos[fcA Tm +ip-{k- 1)A P ] 



3. The Hannay angle in the restricted circular three-body problem 

3.1. Correct derivation of the Hannay angle 

In this section we derive the expression of the Hannay angle, starting from the correct 
representation of the RTBP (see sect. 2.2). Because the RTBP has three degrees of 
freedom, the use of the adiabatic approach is rather cumbersome. Therefore we prefer 
to use a perturbative approach. As shown in the example of the previous section, 
with a perturbation approach we can easily isolate the term that corrects the orbital 
frequency of the small body, at first order in the slow frequency e of P. This term, 
evaluated in the following, is the time derivative of the Hannay angle. 

Although we are ultimately interested in evaluating the Hannay angle for a 
circular motion of T, the eccentricity terms in the perturbation function of the RTBP 
cannot be neglected a priori. Indeed, in section 2.1 we have seen that the Hannay angle 
arises from the quadratic term in P mass, which is generated from the Poisson bracket 
between the original perturbation and the generating function. In the Poisson bracket 
operation, linear terms in the eccentricity (i.e. proportional to Vw) can produce 
eccentricity independent terms. On the other hand, terms of higher order in the 
eccentricity cannot produce such terms. Therefore, it is necessary, and sufficient, to 
consider the expansion of the original Hamiltonian up to order 1 in only, namely: 

H = -tt^t + eA P + m P J2 {h^r) exp[ifc(A T m - A Pm )]+ (30) 

kez 



2A| 



c+(A T )V*exp{t[fcA T m - (k + l)A Pm - tp}} + c k (A T )V* exp{t[fcA Tm ~(k- l)A Pm + 

where Ap is the mean longitude of P, Ap is the conjugate action of P (identical to 
true longitude for a circular orbit), rap is the mass of P relative to A and e is the 
orbital frequency of the P relative to T. In l|3(J|l we have assumed G = Ma — 1 and 
omitted these quantities for simplicity of notation. The symbol l denotes the imaginary 



unit 



T). For d'Alembert rules, the coefficients b 



k,c 



k > ^k 



are real and bu = 6_ 



— c_ k (see ch. 1 of 6 ). In (|30|l the motion of Ap is assumed to be linear in time, 
which reflects the linear dependence of H on Ap. As in the previous section, terms 
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independent of the angles Ap m i Ap, ip are grouped in the integrable approximation Ti. 
and the harmonic terms constitute the perturbation Ji\. 

Using J7J and the new expressions for 7i and 7i\, the function \ required to 
eliminate TL\ from the new Hamiltonian is: 



ib k (A) 



' 2 "I I -p- - e ! k 



• exp[ik(X T - Ap)] + 



A 3 

1V Tr, 



- (fc + l)e 



exp{t[fcATm - (k + l)Ap - 



(31) 



t.c. 



'(A T )Vf 



A3 
L J1 Tm 



-(A-l)e 



exp{i[fcArm - (fc - 1)A P + 



where for simplicity we have omitted the primed notation for the new variables. This 
function is well defined if l/A^ ^ e, and if l/A^ e(k ± 1), which is true away from 
the 1/1 mean motion resonance and from all mean motion resonances of first order in 
the eccentricity. 

The terms of the new Hamiltonian that are of order 2 in Mp can be computed by 
the left hand side of the expression i|l(J|l . After that the Poisson brackets are calculated 
in 110|) , we can finally neglect the resulting terms that depend on the eccentricity, and 
retain only those of order in e (i.e. the terms independent of which are: 



TC'n 



M% v v- - 



2^ 2^ dk,i exp[t(fc + l)(X T „ 
k i^o 



with 



dk.i — — 



db k 



dAi 



A 3 T 



kb k 



db, 



Ap)] 



36, 



dA 7 



A 3 



(32) 



(33) 



-k W 



-k W 



As in the previous section, we now average the new Hamiltonian over Apm- Of 
all the terms in (|32J) only those with k = —I survive, because they are independent of 
Apm- The expansion of these terms at order 1 in e gives the function 



k=ta 



-2b, 



^ + (4)2(fc± l) 



dA T 



k' 2 



k 2 



6A 5 T b 2 k 



(34) 



The change in the frequency of A^ m due to the slow motion of Ap is then 



dA T 

and the Hannay angle becomes 
Oh = — ^Ap m 



(36) 
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3. 2. Correct expression for the Hannay angle 

In order to apply the results of the previous section to the perturbation (J37J, limited 
to the power 2 of Rta I Rap, we set 

fc = L2 = l^ (37) 
for the terms independent of the eccentricity, and 



- 1 4 Ma lot 



8 Af/i lot 



c - = c ^=-l±^ (40) 
8 Af^ lot 

for the terms proportional to the eccentricity. Using (|34[l -il3t) [) we get the following 
contributions to the Hannay angle: 

^) = _^wVm 6 (4i) 

for the terms not related to the eccentricity and 

s<«> = i.m„ f 



"« = 130 'l^J fcj (42) 

for the terms coming from the eccentricity. Note that the major contribution of the 
last term which changes the sign of the effect. In [5], there are two errors concerning 
the sign. Firstly, the sign of the total effect is indeed positive and not negative. 
Secondly, the negative sign would be equivalent to a lengthening of the orbit and not 
to shortening as stated before eq.49 in 0. 

Finally, the total Hannay angle due to P is given by: 



16 " \M A J \Rap J ^ 
The Hannay angle of the Earth due to the perturbation by Jupiter amounts to 
1.14 x 10~ 8 radians which is 1705 m per Jupiter orbit, i.e. 144 m per year, of forward 
displacement, or 4.8 ms of orbital period shortening per year. Furthermore, the value 
of the Hannay angle is comparable (apart the sign) to the value found in even if 
the ratio Rta I Rap is to the sixth power rather than fourth as in since there is A 
large numerical factor (1261/16), eq03l instead of 3 in 

The Hannay effect on an artificial satellite adiabatically driven by the Moon as 
indicated by jS] is more complex, especially if coupled to the issue of measurability, 
and it deserves a separate publication. 



4. Conclusions 



The relation of the Hannay angle to classical perturbation theory of celestial mechanics 
has been given, for both action-dependent and independent perturbations. The 
expression of the Hannay angle in a restricted circular three-body problem has been 
properly derived. 
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Appendix for the three-body problem 

The total potential energy of the system is |S]: 

(M A M P M T M P M T M A \ 

E v = -\ — + — + — G (44) 

\ Map ttTP Mta J 

where R is the interdistance among bodies: 

Rap = [(x A - x P ) 2 + (y A - y P f + (z A - z P ) 2 } 1/2 (45) 

Rtp = [(x T - x P ) 2 + (y T - VP? + (z T - z P ) 2 ] 1/2 (46) 

Rta = [(it - x A ) 2 + (y T - yA) 2 + (zt - z A ) 2 ] 1/2 (47) 

The equations of motion of A and T are, respectively: 

•• dE v w dE v w .. dE v 

M T x T = -—^- M T y T = -—^ Mtzt = ^ — (48) 
oxt oyr ozt 

■■ dE v „ ■■ dE v , r .. dE v tAnS 

M A x A = - — M A y A = -— - M A z A = ~- — (49) 
ox A dy A dz A 

For (^, rj, C) relative coordinates of T and P, relative to A: 

X T = X A + £r Vt — VA + VT ZT = z A + Ct 

xp = x A + S,p yp = yA + r)p z p = z a + (p 

we have: 

1 dE v 1 dE v 



£t - 

Mt dxT M A dxA 
M P (x T -xp) M A (x T -x A ) , Mp(x A -x P ) i M T {x A -x T )] G 



n3 c>3 e>3 r>3 

TP T A AP T A 
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Setting 



we finally have: 



(M T + M A )(, T Mp(6t-6>) M P £p 



n TA 



R AP 



G 



V = 



Mr 



Rtp R\ p 



(6t£p + VtVp + (t(p) 



G 



-{M T + M A )G' 
-{M T + M A )G 
-(Mt + M A )G 





dV 


7? 3 

TA 




rfr 


dV 


Ft 3 

TA 


drjT 


Ct 


l dV 


n TA 





For a angle between T and A: 

1 Rta 



V = GM P 



Rtp R? 



ap 
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(50) 

(51) 

(52) 
(53) 
(54) 

(55) 



and using Legendre polynomials P n (x), the inverse of the distance between P and T 
is written as: 

(56) 



Rtp (ij2 a + R 2 Ap _ 2R TA R AP cos a) 



1/2 



Ra 



1 oo 

-Y, p - 



( Rta \ 



1 



71=0 



V Rap J Rap 



1 , R- TA , 1 to 2 1 \ ( ^ta \ 

1 + — — cos a + -(3 cos a — 1) — — 



Ra 



\Rap J 



Substituting eq.JSSJ) into ca. (|55|l . gives rise to ea. (|T%)l . 



